On-off intermittency over an extended range of control parameter 
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We propose a simple phenomenological model exhibiting on-off intermittency over an extended 
range of control parameter. We find that the distribution of the 'off' periods has as a power- 
law tail with an exponent varying continuously between —1 and —2, at odds with standard on-off 
intermittency which occurs at a specific value of the control parameter, and leads to the exponent 
—3/2. This non-trivial behavior results from the competition between a strong slowing down of the 
dynamics at small values of the observable, and a systematic drift toward large values. 
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Numerous non-equilibrium systems exhibit an irreg- 
ular dynamics, in which quiet periods alternate with 
bursts of activity, giving rise to on-ofF intermittency. This 
phenomenon has been characterized experimentally in 
many different systems, ranging from electronic devices 
spin- wave instabilities plasmas [3], liquid crys- 
tals 0, multistable fiber laser iQ], and nanoscopic 
systems such as nanocrystal quantum-dots [Tj, or dif- 
fusing nanoparticles Q. A standard characterization of 
the intermittency phenomenon is the distribution of time 
spent in the off-state. In most cases, this distribution ex- 
hibits a power-law regime with an exponent close to —1.5, 
with a cut-off at large times [H-Q- From a theoretical 
viewpoint, on-off intermittency is understood as result- 
ing from the effect of multiplicative noise on a system 
close to an instability threshold, generating an effective 
instability threshold that fluctuates over time, and lead- 
ing to the alternation of stable (off, or inactive) and un- 
stable (on, or active) periods. According to this scenario 
[ol-fT^ , the distribution of the duration of the off-periods 
has a power-law tail with a well-defined exponent —3/2. 
This exponent value has a simple interpretation in terms 
of first-return time statistics: on a logarithm scale, the 
dynamics can be mapped onto a random walk on a half 
axis (the one corresponding to small values of the physical 
variable) , so that the duration of off-periods corresponds 
to the first return time to the origin of a random walk, 
whose statistics is known to behave with a tail exponent 
-3/2 13. 

In spite of the large body of results corroborating this 
scenario, some experimental systems have been found to 
exhibit a different behavior, with exponent values be- 
tween — 1 and —2 In addition, it has also been 
observed experimentally in some cases that the power-law 
distributions of the off-periods appear for an extended 
range of control parameter |l^, and not only close to 
a given threshold. While a possible explanation of the 
variations of the exponent might be the presence of cor- 
relations in the dynamics jTZj, |1_S] , no theoretical scenario 
is, to our knowledge, presently able to account for the 
existence of on-off intermittency over a finite range of 
control parameter. 

In this note, we study both analytically and numer- 
ically a simple stochastic model exhibiting on-off inter- 



mittency, with a continuously varying tail exponent, over 
an extended range of control parameter. We also provide 
a simple interpretation for the upper and lower bounds 
of this range, and give an analytical argument, based on 
random walk properties, to determine the tail exponent 
as a function of the control parameter. Intermittency re- 
sults in this case from a competition between a slowing 
down at small values of the observable, and an average 
drift toward larger values. 

Model and dynamics - We consider a simple model de- 
scribing the stochastic dynamics of a continuous positive 
variable x, corresponding to an arbitrary physical observ- 
able. The dynamics is defined as follows. Starting from 
a value x, the process can jump to a value in the interval 
[x' , x' + dx'] with a probability per unit time iy(x'|a:) dx' , 
where W^(a;'|a;) is given by 



Wix'\x)^iyopix')ij{ - 
' x 



(1) 



vq being a characteristic frequency. The distribution 
p{x'), normalized according to p{x') dx' = 1, charac- 
terizes the a priori probability to choose the value x' (it 
may thus be thought of as a density of states), while 'ip{u) 
describes an acceptance rate of the transition from x to 
x' . This rate is chosen to be a function of the ratio x' /x, 
so that the process can be considered as multiplicative. 
The density p{x) is assumed to behave for small values 
of a; as a power law, namely 
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where 70 and c are positive constants. For reasons that 
will appear clear later on, the function ^{u) is assumed 
to fulfill the symmetry 



ili{u) 



u> 0, 



(3) 



with 7 > a parameter of the model. When an explicit 
form of is required, we use in the following 
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FIG. 1: Examples of trajectories for increasing values of 7. 
(a) 7 = 0.5, (b) and (c) 7 = 1.5, (d) 7 = 2.5 (ry = 0,vo = l). 
(c) corresponds to the same data as (b), but plotted on a 
logarithmic time scale. 

^ 0)j which fulfills the symmetry property Eq. ([3]). In 
what follows, we consider 70 and rj as fixed parameters, 
and 7 as the control parameter of the system. 

On- off intermittency: qualitative approach.- We show 
on Fig. [1] examples of trajectories corresponding to dif- 
ferent values of 7, using the function ipex(u) defined in 
Eq. (HI), with 77 = 0, and p{x) = 70^'^''"^ (0 < a; < 1). We 
see on Fig.[TJa) that for small values of 7, the process re- 
mains in a fluctuating (or active) phase. In contrast, for 
large enough values of 7, the dynamics quickly converges 
to X = 0, as seen on Fig.[IJd). For intermediate values of 
7, on-off intermittency appears, with a succession of ac- 
tive periods separated by quiet periods of time in which 
X remains very small [Fig. [Ijb)]. On a logarithmic time 
scale [Fig. [IJc)], one observes that lower and lower val- 
ues of X are progressively reached, while bursts of activity 
where x ~ 1 remain present even at large time. 

The fact that the process remains trapped for long pe- 
riods of time at small values of x can be understood from 
the transition rates The average sojourn time Ts{x) 
at a given value x can be evaluated by integrating over 
all possible escape paths: 



1 



Ts{x) 



dx'W{x'\x). 



(5) 



From the small x expansion of p(x) given in Eq. ([2]), 
one finds Ts{x) ^ tix~'^° for x — > 0, with ti — 
[vqc dv v'^°^^^j{v)\^^ , leading to a divergence of Ts{x) 
in this limit. 

As a first characterization of intermittency, we deter- 
mine the stationary distribution pst{x). Using the sym- 
metry ([3]), one finds that a detailed balance relation [l^ 
W {x' \x) Pst{x) = W{x\x')pst{x') holds with respect to 
the distribution 



(6) 



implying that Pst{x) is the stationary distribution. We 
however need to check that the normalization constant 



Z = Jp°° dx p(x) is finite, otherwise the distribu- 
tion Pst{x) does not exist (it cannot be normalized). 
The convergence at the upper bound is ensured, but 
the integral may not converge at its lower bound. Us- 
ing the small x asymptotic behavior of p(x), we find 



p{x) X 



cx 



70-7-1 



so that the integral defining Z con- 



verges only if 7o — 7 > 0. Hence the distribution Pst{x) is 
well-defined for 7 < 70 but it becomes non-normalizable 
for 7 > 7o, which indicates that the steady-state distri- 
bution should in this case be a Dirac delta function at 
a; = For 7 > 70, the time spent around a; = thus 
completely dominates the dynamics, which becomes non- 
stationary since the time needed to reach the asymptotic 
delta distribution is infinite. 

To sum up, the fact that Pstix) becomes non- 
normalizable accounts for the change of behavior be- 
tween Fig. [TJa) and (b)-(c). We now need to understand 
the mechanism responsible for the different behaviors ob- 
served in Fig. mb)-(c) and (d). In order to identify the 
value of 7 separating these two regimes, we determine 
the average trend of the process on a logarithmic scale, 
namely (Alna:)^; with Alnx = Inx' — Inx. The average 
is computed over x' for a fixed value of x: 

/>oo 

{A\nx)^= dx' (j){x'\x){\nx' -\nx) (7) 
Jo 

where (l){x'\x) is the probability to jump from x to x' , 
obtained by normalizing the transition rate W{x'\x): 
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From the small x behavior ([2]) of p{x), we get for small 
values of x and x' 



4,ix'\x) = -}(- 
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with the scaling function (f>{u) ~ vqcti u''" ^ip{u). Using 
the symmetry of ipi''^) given in Eq. ([3]), one finds 



{Alnx)x = VoCTi 



— Inw V'(^) i'v'^" 

V 



(10) 



We assume that the integral in Eq. (ITU)) converges, which 
is true in the example of ipcx{u) as soon as 7 > 70 — 77. 
We first note that, in the considered small x limit, the 
expression of (Alna;)^. becomes independent of x. Sec- 
ond, one readily sees from Eq. ([T0|) that {A\nx)x is 
positive for 70 < 7 < 270, and negative for 7 > 279. 
Hence, when 7 > 270, every step tends, on average, to 
reduce the value of x, leading to a roughly monotonous 
convergence of x towards [Fig. (Hd)]. In contrast, for 
7o < 7 < 270, every step increases, on average, the value 
of X, which favors finite values of x instead of small ones. 
Yet, rare negative steps play an important role due to 
the strong slowing down of the dynamics for small values 
of X [Fig. [Ijc)]. The competition between positive drift 
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FIG. 2: Distribution P(t) of the durations r of the off-periods 
for different values of 7. The distribution has been multiphed 
by T to enhance readabiUty. From bottom to top: 7/70 = 1-1 
to 1.9 by steps of 0.1 (7; = 0, vq = 1). 



and trapping at low values leads to the observed on-off 
intermittency phenomenon, which is thus expected to ex- 
ist over a finite range of the control parameter 7, namely 
7o < 7 < 270. Interestingly, the above analysis pro- 
vides an interpretation of the upper and lower bounds, 
7min = 7o and 7max = '^lo, of the range over which on- 
off intermittency is present: above 7inin, kinetic trapping 
effects are strong close to 2; = (the distribution Pst{x) 
becomes non-normalizable), while below 7max, the evolu- 
tion is biased toward large values of x. The existence of 
an overlap between these two ranges results in the pres- 
ence of on-off intermittency. Note that 7niax plays a role 
similar to that of the instability threshold in the standard 
on-off intermittency scenario. 

Distribution of off-periods - In order to characterize 
more quantitatively the intermittent regime, we have de- 
termined numerically the distribution P{t) of the dura- 
tions T of the off-periods, by measuring the first return 
time to a given threshold xth (we have checked that the 
results do not significantly depend on the precise thresh- 
old value as long as Xth ~ 1). Simulations were performed 
using p{x) = 700;'''''"^ and the function ipcxiu) defined in 
Eq. (HI), for different values of rj. The distribution P{t) 
is plotted on Fig. [2] for different values of 7 in the range 
7o < 7 < 270. A power-law behavior is observed in the 
tail of the distribution, namely 



P(r) 



(11) 



for large r. The exponent a is seen to be independent of 
rj, and varies almost linearly with the control parameter 
7 (see Fig. [3]). Small deviations from the power law can 
be seen for 7 close to 270, but the observed variation of 
the exponent a when varying the measurement window 
over a reasonable range remains small, at most of the 
order of the symbol size on Fig. [3] 

We now give an analytical prediction of the exponent 
a. As the sojourn time Ts{x) at small values of x becomes 
very large, it is natural to assume (at variance with usual 



FIG. 3: Exponent a of the tail of the distribution of the off- 
period durations as a function of the ratio 7/70, for different 
values of ry. The straight line is the prediction a = 2 — 7/70. 



first return problems of random walks) that the return 
time r is dominated by the largest sojourn time in the 
trajectory, that is, the sojourn time Ts{xra) at the small- 
est value of Xm reached. Using the small x expansion 
of Ts{x), we thus have the scaling relation t ^ tix:^-^'^°. 
Here again, it is convenient to work with logarithmic vari- 
ables. We thus define the variable z = — 7olna;, from 
which the relation r ^ ri e^'" follows, where Zm is the 
maximal value of z along a trajectory between two cross- 
ings of the threshold Zth — ~7oln2:th- The distribu- 
tion P{t) can be deduced from the distribution ^(zm) 
through P{t) = P{z^)dz^/dT. From the scaling rela- 
tion r ~ Tl e^" , one obtains 



1 ~ / T 
P{t) w -P In — 

T \ Tl 



(12) 



In terms of number of steps (instead of real time), the 
trajectory of the logarithmic variable z is a random walk 
with a step statistics $(?/) obtained from (j){x'\x) in the 
small X limit [see Eq. ([5])] through the correspondence 

y = -70 ln(2:7a::): 



he 
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y > 
y<0 



(13) 



where ^ = 1-f 77/70, /? = (7+?7)/7o-l and 6 = (3 ^j. / {[5 + ^j.) . 
We assume that /i > /?, or cquivalcntly 7 < 270. The 
distribution P(zni) can be determined from the auxiliary 
problem of the random walk (|13|) with absorbing bound- 
aries at z = and z = z^. We introduce the probability 
(5(zth, Zm) that the random walk, starting from a posi- 
tion zth (0 < zth < Zm), is absorbed at the boundary 
Zin. The distribution P{zm) can be determined by not- 
ing that the walks having a maximal value between Zm 
and Zm + dztn are precisely the walks that are absorbed 
by a boundary at Zm, but not by a boundary at Zm + dzm, 
in the auxiliary problem. It follows that 



P(^m) 
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The probabihty Q{z, z,„) satisfies the integral equation 



Q{z,Zm) 



dy^iy)Qiz + y,z„ 



(15) 

Looking for the solution of Eq. (|15p as a sum of exponen- 
tials of z, one finds 



Q{z., Zm) 



(^^e'^P'-^)- - py (16) 



Using Eq. ([Tl)) . we get for large z^ 

P(z„,)«i^(zth)e-(''-'^)^" 



(17) 



with i^(zth) = - /3)[Aie(^~'^)^"' - /3]//x2. From 
Eq. p^ . we obtain for P(t) the expression: 



P(r) 



i^(2th)r[ 



yielding a ~ fj, — f3, or in terms of 7, 



00, 



J7 
70' 



(18) 



(19) 



This value, which is independent of 77, is in good agree- 
ment with the numerical results shown on Fig. [31 

The range 70 < 7 < 270 over which intermittency ap- 
pears corresponds to < a < 1, that is to a distribution 
P(r) with an infinite mean value (r), as in standard on- 
off intermittency for which a = 1/2. If 7 < 70, the 
average duration of off-periods is finite, while if 7 > 270, 
the process becomes completely trapped in the off state, 
without any return to the threshold value, so that the 
distribution P(t) can no longer be defined. We also note 
that the exponent a is independent of the threshold value 



Zth [which only appears in the prefactor K{zth) 
observed numerically. 



as also 



Discussion.- In this note, we have considered a simple 
model exhibiting a finite range of control parameter over 
which on-off intermittency is present, with a continuously 
varying exponent a characterizing the distribution of the 
duration of off-periods. The mechanism at the origin of 
the intermittency phenomenon and of the power-law dis- 
tribution P(t) differs from that found in standard on-off 
intermittency. In the standard case, the statistics of off- 
periods is related to the distribution of return times of 
unbiased homogeneous random walks (whence the expo- 
nent a = 1/2 follows), while in the present model, on-off 
intermittency results from the competition of bias and 
slowing down effects, leading to a non-trivial statistics 
with a tunable exponent a, and to an extended range of 
intermittency. 

Further work is needed to identify more clearly the re- 
lation between the present model and the standard sce- 
nario of on-off intermittency. For instance, it would be of 
interest to make contact with dynamical systems theory, 
for instance by finding models having a behavior simi- 
lar to the present one, but being defined by a stochastic 
differential equation or by a chaotic map. Another open 
question is the physical interpretation of the control pa- 
rameter 7, which is a phenomenological parameter at this 
stage. As a first attempt in this direction, we note that 
the present model can be mapped, for ry = 0, onto the 
Barrat-Mezard model [l^, [2l| , a simple model exhibiting 
aging dynamics. In this mapping, 7 is found to be pro- 
portional to the inverse temperature and thus acquires a 
simple physical interpretation. It would be interesting to 
find other examples of physical realizations of the present 
model. 
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